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The optoelectronic properties of single- and double-wall carbon nanotubes (CNTs) noncovalently 
functionalized with tetraphenylporphyrins (TPPs) are addressed by dispersion-corrected ab initio 
calculations. Five CNT species with different chiralities were considered. We find that the most 
stable configurations are those where the CNTs are fully covered by TPPs, exhibiting binding energy 
of about 2 eV/TPP. The semiconducting CNT-TPP compounds show optical response characterized 
by a strong absorption associated to the TPP bands, with increasing intensity with the TPP con¬ 
centration. In addition, molecular dynamic simulations show that the compounds would be stable 
at temperatures as high as 100°C. 


I. INTRODUCTION 

The noncovalent functionalization of carbon nanotubes 
(CNTs) with porphyrin chromophores has proven to be 
an efficient method that combine the strong photoab¬ 
sorption of porphyrin with the CNT electron transport 
and charge separation properties. This supramolecular 
assembly has been explored recently as light-harvesting 
donor/acceptor complexesi^^— , as well as for the electro¬ 
catalysis of the oxygen reduction reactionii^— . Some 
schemes proposed for the such a functionalization in¬ 
clude the CNT encapsulation in micelles of surfactanU^ 
or in polymers. Recently, CNT-porphyrin complexes ob¬ 
tained from a micelle-swelling technique have shown an 
extremely efficient energy transfer from the photosensi¬ 
tized porphyrin to the nanotubes with a remarkable sta¬ 
bility, confirming the potentiality of these complexes for 
light-harvesting applications^^. However, some key pa¬ 
rameters like the TPP concentration on the CNT surface 
and the stability of the compounds are not yet clarify. 

Another important issue concerning the noncovalent 
functionalization of CNTs with organic molecules or 
polymers is the ability to discriminate between nan¬ 
otube species in terms of either diameter or chiral angles. 
Experiments have reported a substantial separation of 
single-wall CNT according to their chirality (metallic ver¬ 
sus semiconducting) after the adsorption of porphyrins 
and organic polymer o^^d^ . Other experiments reported 
improvement in the preparation of CNT-polymer solu¬ 
tions with a high degree of selectivity, favoring CNTs 
with large chiral angles and the inhibition of the metal¬ 
lic species Although the mechanism causing the selec¬ 
tivity is not clearly understood, the above experiments 
suggest that the CNT-polymer bonding is strongly in¬ 
fluenced by both the relative orientation of the poly¬ 
mer chain to the nanotube structure, as well as possible 
charge transfer from metallic CNTs. 

In the present work we investigate semiconduct¬ 
ing single- and double-wall CNTs noncovalently func¬ 
tionalized with free-base tetraphenylporphyrin (TPP) 
molecules. Our goal is to give insight in the stability 
of such complexes with temperature and also determine 


its electronic properties and optical response. Addition¬ 
ally, we consider both metallic and semiconducting CNTs 
with a similar diameter in order to identify a possible chi¬ 
rality dependence in the TPP adsorption strength. Our 
results show that the TPP molecules tends to cover the 
whole CNT surface, driven by quite strong dispersive 
forces. The CNT-TPP binding energies would depends 
on the CNT diameters instead of chirality, as previously 
found for a single TPP adsorbed on different CNTs^S. 
We also find that the CNT-TPP complexes would be 
stable at temperatures as high as 100°C according to 
our molecular dynamic simulations. Interestingly, the 
TPP absorption bands remain almost unchanged after 
the TPP stacking on the CNT surface, showing a small 
redshift. However, the intensities of the TPP absorption 
bands increase with the concentration of molecules on the 
CNTs, independently if they have single or double wall. 


II. COMPUTATIONAL METHOD 

The calculations are carried out using density func¬ 
tional theory (DPT) including Van der Waals (VdW) 
dispersion-corrected exchange-correlation functional, as 
implemented in the SIESTA ab initio package^i. A 
double-C, singly polarized basis of localized atomic or¬ 
bitals is used, while the core-electron interaction is de¬ 
scribed with norm-conserving pseudopotentials. We in¬ 
vestigate the TPP adsorption on the metallic (10,10) and 
semiconducting (16,0) single-wall CNTs, with diameters 
of 13.9 and 12.8 A, respectively. We also consider the 
double-wall CNTs (10,10)@(5,5) and (16,0)@(8,8), where 
the inner CNTs have diameters of 7.0 and 6.4 A, respec¬ 
tively. These CNTs can host up to five TPP molecules 
around the perimeter. In order to describe the CNTs 
fully covered with a single layer of TPPs, we use a su¬ 
percell with a volume of (38 x 38 x oq) A^, and periodic 
boundary conditions along the nanotube axis {z direc¬ 
tion). For (16,0) and (10,10) CNTs, oq is chosen to be 
four and seven times de CNT lattice constant, of 4.334 
and 2.500 A, respectively. For the double-wall CNTs, 
the oo parameter is the same. The cell size along z can 
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host one molecule, giving a separation between images 
of about 3.8 A. These supercells allow us to describe the 
(10,10)-5TPP and (16,0)-5TPP compounds with a max¬ 
imum concentration of about 0.3 TPP/A. 

We also study the (6,5) CNT fully covered with a sin¬ 
gle layer of TPPs. This CNT has a diameter of 7.7 A 
and can host up to four molecules around the perime¬ 
ter, as previously established^. To describe the CNT 
periodicity along z, we use a supercell with a volume of 
(30 X 30 X ao) A^, where ao is the CNT lattice constant, 
of 41.327 A. The cell size along z can host a maximum of 
two molecules side by side, therefore with this supercell 
we can describe two adjacent rings of TPPs around the 
CNT, that is the (6,5)-8TPP compound. This configu¬ 
ration represents the full surface coverage with a concen¬ 
tration of about 0.2 TPP/A. 

The self-consistency of the density matrix is achieved 
with a tolerance of 10“"^ and a common basis size or en¬ 
ergy shift of 10 meV is applied. A grid cutoff of 100 Ry 
and the P point were used for the real-space and k- 
space integration, respectively. The accuracy of these 
parameters were tested considering larger grid cutoff of 
150 Ry and k mesh of 1 x 1 x 10. Negligible variation 
in total energy, optical spectra and band structures were 
found, ensuring converged results. For the optical cal¬ 
culations, we use a 1 X 1 X 31 A:-point mesh, for inci¬ 
dent light polarized along the CNT axis. The CNT-TPP 
binding energy is calculated by the energy difference be¬ 
tween adsorbed and separated constituents, considering 
corrections due to the basis set superposition error. The 
CNT-TPP systems were fully relaxed by conjugate gradi¬ 
ent minimization until the forces on the atoms were less 
than 0.05 eV/A. The dynamical properties and stability 
of the (16,0)-TPP, (10,10)-TPP and (6,5)-TPP systems 
were investigated by molecular dynamic simulations in 
the NVT ensemble using the Nose-thermostat approach 
at 373 K over a total simulation time of 4 ps, with a time 
step of 1 fs. 

The physisorption of the TPP molecules on the CNT 
surfaces is assessed by VdW density functional. We test 
different parametrizations for the VdW functional by cal¬ 
culating the benzene adsorption on graphene. These 
parametrizations include those of Dion et al. (DRSLL)^; 
Klimes et al. (KBM)^; Vydrov and Van Voorhis (VV)^; 
Cooper (C09)M; Berland and Hyldgaard (BH)^; and Lee 
et al. (LMKLL)^. Figure [T] shows our results for the in¬ 
teraction energy as a function of the adsorption distance 
for the different VdW parametrizations. The square box 
indicates the experimental binding energy of benzene on 
graphite of 0.5 ± 0.08 eV— , whereas the variation in dis¬ 
tance ranges between the experimental interlayer sep¬ 
aration of graphite (3.35 A) and that of double layer 
graphitic carbon (3.84 A)^. As can be seen in Fig. [H 
the parametrization that best approaches to the experi¬ 
mental binding energy of benzene on graphite is that of 
DRSLL. The difference with respect to the experimental 
results is about 0.15 eV, which in some way represent 
the error of the binding energy calculations. The VdW 


parametrization of DRSLL will be used throughout this 
work. 



FIG. 1. Interaction energy as a function of the separation for 
the benzene adsorption on graphene calculated with different 
parametrizations for the Van der Waals density functional. 
The square box indicates the available experimental data. 

The optical response of the functionalized CNTs is ob¬ 
tained through first-order time-dependent perturbation 
theory, by calculating the imaginary part of the dielec¬ 
tric function (£ 2 )- £2 gives us a first approach for the 
optical absorption coefficient and it is calculated accord¬ 
ing to the equation: 

£ 2 ( 0 ;) = A f dk^|e^(^,(k)|r|^„(k))|2 5(^;e-^;„-M, 

C,V 

( 1 ) 

where A is a constant that depends on the cell size; 
and 'I'^ are the occupied and empty Kohn-Sham orbitals, 
respectively. The delta function represents the conserva¬ 
tion of energy, which is described by a gaussian function 
with a smearing of 0.06 eV. We must note that DFT 
calculations do not include electron-hole and electron- 
electron interactions, the so called many-body effects. 
An accurate description of the optical response or pho¬ 
toexcitations in CNT-TPP compounds needs theories be¬ 
yond DFT. However, these calculations are currently pro¬ 
hibitive considering the size of the systems under study. 
Thus, £2 as calculated by perturbation theory can give 
a first approach for the CNT-TPP absorption spectra, 
which are typically redshifted by about 0.8 eV with re¬ 
spect to those including many body effects. 

III. RESULTS AND DISCUSSIONS 
A. Stability and electronic properties 

Figure [H illustrates the TPP molecules adsorbed on 
the semiconducting (16,0) CNT in the equilibrium geom¬ 
etry. The molecules are initially placed in positions close 
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to those depicted in Fig. [2l but with larger separations 
among them. Then, the whole system is allowed to relax 
without any constraint. The most stable configuration 
is found for an array of five TPP molecules covering the 
CNT perimeter, forming a ring around the CNT. The 
system evolves in such a way that the TPP molecules 
tend to fit in the most compact configurations, that is 
with the phenyl group of one molecule between those of 
the neighboring molecules. The TPP binding energy in 
the (16,0)-5TPP complex is found to be of 2.10 eV/TPP, 
which is 0.14 eV higher in energy than the adsorption of 
a single TPP, suggesting that the full surface coverage 
is energetically favorable. The TPP adsorption on the 
double-wall CNT (16,0)@(8,0) shows a binding energy of 
2.20 eV/TPP. This suggests that the relevant interaction 
behind the strong adsorption occurs between the adja¬ 
cent aromatic systems with a small contribution from the 
inner CNT. For the metallic compounds (10,10)-5TPP 
and (10,10)@(5,5)-5TPP, we observe a similar behavior. 


(a) (b) 



FIG. 2. Schematic representation of TPP molecules adsorbed 
on single- and double-wall CNTs in the equilibrium geometry. 
(a)-(b) Top and front views of (16,0)-1TPP, (c) (16,0)-5TPP, 
and (d) (16,0)@(8,0)-5TPP. 

Previously, we reported calculations for the TPP ad¬ 
sorption on the semiconducting (6,5) CNT, of 7.7 A in di¬ 
ameter, considering different concentrations of adsorbed 
molecules^^. The most stable compound is found for a 
configuration of four TPPs covering the CNT perimeter, 
with a binding energy of 2.45 eV/TPP. In the present 
work we have realized that the TPP binding energy 
strongly depends on the parametrization of the VdW in¬ 
teraction, as shown in the Fig. [TJ Thus, we have reeval¬ 
uated the TPP binding energy for the (6,5)-4TPP com¬ 
pound, considering the best parametrization for the VdW 
density functional. Our results show a lower binding en¬ 
ergy of 1.79 eV/TPP, which we believe is more realistic 
than our previous result. 


Table U lists the adsorption distances and binding en¬ 
ergies for the TPP molecules adsorbed on the single- and 
double-wall CNTs. According to this results we can infer 


TABLE I. Adsorption distance (r) and binding energy (Ei,) 
of TPP molecules adsorbed on single- and double-wall CNTs. 
D (d) is the diameter of the outer (inner) double-wall CNTs. 


Compound 

D (A) 

d(A) 

r (A) Eb (eV/TPP) 

(6,5)-lTPP 

7.68 

— 

3.64 

1.31 

(6,5)-4TPP 

7.65 

— 

3.14 

1.79 

(6,5)-8TPP 

7.65 

— 

3.14 

1.95 

(16,0)-1TPP 

12.81 

— 

3.16 

1.96 

(16,0)-5TPP 

12.82 

— 

3.23 

2.10 

(16,0)@(8,0)-1TPP 

12.87 

6.44 

3.12 

1.97 

(16,0)@(8,0)-5TPP 

12.92 

6.44 

3.18 

2.20 

(10,10)-1TPP 

13.89 

— 

3.24 

2.02 

(10,10)-5TPP 

13.88 

- 

3.23 

2.26 

(10,10)@(5,5)-1TPP 

13.90 

7.02 

3.20 

2.04 

(10,10)@(5,5)-5TPP 

13.92 

7.02 

3.21 

2.24 


some general properties for the CNT-TPP compounds: 
(i) The TPP binding energy tends to increase with the 
CNT diameter, somewhat expected due to the CNT- 
TPP TT-stacking interaction, which is maximum for pla¬ 
nar systems, (ii) For single- and double-wall CNTs, the 
TPP binding energy is almost the same, suggesting that 
the adsorption interaction is localized, that is, depending 
mostly on the adjacent aromatic systems, (iii) The num¬ 
ber of molecules on the CNT surface tends to increase 
the adsorption energy, favoring the full surface coverage. 
We observe an increase in the compound binding energy 
of about 0.2 eV per molecule for higher concentrations. 
We attribute this result to an additional VdW attractive 
force among the TPP phenyl groups which would play 
a key role in the stabilization of the TPP layer covering 
the CNT. 

Figure |3] shows the density of states (DOS) close to the 
Fermi level for the single- and double-wall CNTs fully 
covered with TPP molecules. The Fermi energy {Ep) is 
placed at the center of the energy gap for the semicon¬ 
ducting systems. The vertical lines in the figure indicate 
the frontier orbitals of the isolated molecule: HOMO-1, 
HOMO, LUMO, and LUMO-l-1. The states added to the 
CNT DOS by the HOMO-1 and HOMO levels are iden¬ 
tify by peaks at around -1.0 and -0.5 eV, respectively. 
Whereas the states associated to LUMO and LUMO-l-1, 
which are very close in energy, appear as a larger single 
peak at around 1 eV. The height of these peaks are re¬ 
lated to the number of TPP molecules adsorbed on the 
CNTs. On the (6,5) CNT we put eight TPPs per super¬ 
cell, referred as (6,5)-8TPP in the Fig. [Sja). Therefore, 
the molecular peaks are larger than those found in the 
(16,0)-5TPP and (10,10)-5TPP compounds, containing 
five molecules, as shown Fig. [3Kb) and 3(c), respectively. 
We verify that the strong CNT-TPP interaction induces 
distortions in the adsorbed molecules, resulting in small 
changes in the position of the energy levels, which ex¬ 
plain to some extent the width of the absorption peaks 
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and also the decrease in the HOMO-LUMO energy with 
respect to the free molecule, as shown in the Fig. [31 





FIG. 3. Density of states for the CNT-TPP compounds: (a) 
(6,5)-8TPP, (b) (16,0)-5TPP and (16,0)@(8,0)-5TPP, and (c) 
(10,10)-5TPP and (10,10)@(5,5)-5TPP. The vertical lines rep¬ 
resent the frontier orbitals of the free TPP molecule. 

For the (6,5) CNT we observe a reduction in its 
band gap as induced by the TPP adsorption, as shown 
Fig.lSKa). This behavior can be explained for the position 
of the TPP HOMO levels which coincide with the CNT 
valence band maximum, inducing the rising of these lev¬ 
els inside the CNT band gap, as previously discussed for 
the (6,5)-4TPP compound^^. Regarding the double-wall 
CNTs fully covered with TPP molecules, we observe es¬ 
sentially the same electronic properties than those of the 
single-wall CNT, unless a small shift to lower energies 
which is more pronounced for the semiconducting CNT 
species, as shown Figs. 3 (b) and 3(c). We attribute 
this shift to a weaker interaction between the molecules 
and the outer CNT due to a screening effect, which are 
also expressed in the larger TPP adsorption distance, as 
shown in Table [H 

B. Optical absorption properties 

The optical absorption spectroscopy of the TPP 
molecule shows resonances at 2.831 and 2.357 eV, which 
are associated to the B (Soret) and the first Q bands, 
respectively!^. B and Q bands split in two electronic 
components, x and y, owing to the symmetry break¬ 


ing in the TPP molecule due to the alignment of the 
N-H bonds in one direction. Theoretically, the opti¬ 
cal response of the CNT-TPP compounds can be esti¬ 
mated through the imaginary part of the dielectric func¬ 
tion (£ 2 ), which is proportional to the optical absorption, 
considering a light polarization parallel to the tube axis. 
For the free TPP, our results show clear resonances at 
2.04 and 1.60 eV, which we associate to the B and Q 
bands, respectively. These resonances can be character¬ 
ized by transitions between TPP frontier orbitals^ Con¬ 
sequently, Qx (Qy) can be associated to the transition 
from HOMO to LUMO (LUMO-l-1), while Bx (By) to the 
transitions from HOMO-1 to LUMO (LUMO-l-1). The 
close proximity in energy between LUMO and LUMO-f 1 
levels, of less than 0.1 eV, leads to the appearance of sin¬ 
gle peaks in the optical spectrum of the compound, which 
we called Qxy and Bxy We observe that the theoretical 
energy difference between Qxy and Bxy bands (0.44 eV) 
is close to those observed in experiments (0.47 eV)!^. 
Therefore we can say that the difference in energy be¬ 
tween Qxy and Bxy is well represented by the frontier or¬ 
bital approximation^. However, experiments show the B 
band at 2.831 eV, while in our theoretical approach it is 
found at 2.04 eV. Thus, our single-particle approximation 
shows a general redshift of 0.8 eV in the absorption spec¬ 
tra of the CNT-TPP compound, but qualitatively pre¬ 
serving the characteristic of the TPP absorption bands. 

It is important to note that excitonic effects alter dra¬ 
matically the optical spectra in semiconducting CNTs. 
For each peak obtained with the single-particle approxi¬ 
mation, it is found a series of visible exciton lines when 
hole-electron interaction is included^. Optical excita¬ 
tion of CNT in photoluminescence or Raman experi¬ 
ments have confirmed these excitonic states, which for 
CNTs with diameters around 8 A lie several hundred of 
meV below the conduction band edge^. If both electron- 
electron and electron-hole interactions were included in 
our calculations, no redshift on the CNT optical transi¬ 
tion should be expected. In addition, because the molec¬ 
ular absorption bands are located far from the CNT 
conduction-band edge, the optical transitions of the ad¬ 
sorbed molecules should not be altered. This has been 
confirmed by optical absorption spectroscopy of isolated 
TPPs and the CNT-TPP compounds^. 

Figure |3| shows the imaginary part of the dielec¬ 
tric function for the semiconducting CNTs functional¬ 
ized with TPP molecules with different concentrations, 
as compared with the pristine CNTs. For the fully 
coverage configurations: (6,5)-8TPP, (16,0)-5TPP and 
(16,0)@(8,0)-5TPP, we clearly identified two peaks at 
around 1.5 and 2.0 eV which are associated to the Qxy 
and Bxy bands of the adsorbed molecules, respectively. 
The intensity of these peaks can be related with the con¬ 
centration of adsorbed molecules. For the (6,5)-8TPP 
compound, we observe a redshift of the Bxy band by 
about 0.1 eV. This redshift can be associated to the 
screening of the Coulomb interactions in the nanotube by 
the adsorbed molecules, in agreement with experimental 
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FIG. 4. Imaginary part of the dielectric function for the semi¬ 
conducting CNT-TPP compounds: (a) The (6,5)-TPP, (b) 
The (16,0)-TPP, and (c) The (16,0)@(8,0)-TPP. The dashed 
lines indicate the Qxy and Bxy absorption bands of the free 
TPP molecule. 

results^. However, the opposite situation occurs for the 
(16,0)-5TPP and (16,0)@(8,0)-5TPP compounds where 
the Qxy band is redshifted while the Soret band remains 
almost unchanged. Additionally, we note that the TPP 
adsorption on the different CNT-TPP compounds appear 
to be not affected by the nature of the CNT, suggest¬ 
ing that multi-wall CNTs would exhibit the same optical 
properties that those observed in single-wall CNTs. 


C. Molecular dynamic simulations 

The stability of the CNT-TPP compounds is investi¬ 
gated through ab initio molecular dynamic simulations at 
100°C, which is a typical temperature for a device opera¬ 
tion. We consider a total simulations time of 4 ps. After 
a stabilization during 2 ps, we analyze the adsorption dis¬ 
tance r of each TPP molecule on the CNT surface for the 
compounds (10,10)-5TPP, (16,0)-5TPP, and (6,5)-4TPP. 
The evolution of the adsorption distance with time is ob¬ 
tained by taking the minimum distance from a TPP N- 
atom to the closer nanotube C-atom, with the help of the 
LPMD code toolkit^. Our results are shown in Fig. [5] 
As can be seen, the TPP molecules fluctuate onto the 
CNT surface without desorbing. The mean adsorption 
distances of the five TPP molecules on the (10,10) and 


(16,0) CNTs are found to be of 3.619 and 3.623 A, respec¬ 
tively, which are 11.6% and 12.2% larger than those found 
for the equilibrium geometry at 0 K. The high TPP bind¬ 
ing strength is also observed in the smaller (6,5) CNT, 
with a diameter of 7.7 A, as shown our MD simulation 
for (6,5)-4TPP. Here the mean adsorption distance, of 
3.550 A, is 13.1% larger than the adsorption distance at 
0 K. This confirms the strong attractive interaction of 
the (6,5)-TPP compounds as experimentally reported^. 
The above results suggest that the TPP molecules remain 
adhered to the CNT surface at temperature as high as 
100°C. However, the large TPP binding energy, of about 
2 eV, suggests a much higher desorption temperature. To 
check this, we performed additional MD simulations for 
(16,0)-5TPP, considering a temperature of 500°C. Our 
results reveal that the TPP molecules still remain ad¬ 
hered on the CNT surface, confirming the high stability 
of the CNT-TPP compounds. 





FIG. 5. Adsorption distance r as a function of time for the 
molecular dynamic simulation of the CNT-TPP compounds 
at 373 K. The dashed lines indicate the average adsorption 
distance. The different curves indicate the evolution of r for 
each TPP molecules. 

We also note that no sizable changes in the adsorp¬ 
tion distance between the metallic and semiconducting 
CNTs is obtained. Considering that the metallic (10,10) 
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and semiconducting (16,0) CNTs have similar diameters 
and they can host up to five TPP molecules around its 
perimeter with almost the same binding energies (see 
Table |I|, our results suggests that the CNT chirality 
would not be a relevant parameter for the TPP binding 
strength. 

IV. SUMMARY AND CONCLUSIONS 

In summary, using dispersion-corrected DPT calcula¬ 
tions and linear optical response we have investigated 
the stability, electronic and optical properties of single- 
and double-wall CNTs noncovalently functionalized with 
TPP molecules. CNTs with different chiralities were also 
considered. Our results show that the most stable con¬ 
figuration for all CNTs under study are those where they 
are fully covered by TPP molecules. This configuration 
shows binding energy of around 2 eV/TPP for all the 
CNT species. We find that the molecular aggregation on 
the CNT surface is an exothermic process with a gain in 
energy of about 0.2 eV. We also find that the complex 
stability would be independent of both the CNT chirality 
and the CNT number of walls. Therefore, no evidence of 
a preferential interaction of TPP molecules with metal¬ 
lic or semiconducting CNTs is found. However the CNT 
diameter appears to be the relevant parameter for the 
increasing TPP binding strength. 

The semiconducting CNT-TPP compounds exhibit op¬ 


tical response characterized by a strong adsorption asso¬ 
ciated to the TPP bands, with increasing intensity with 
the number of adsorbed molecules. This characteristic 
is preserved in double-wall CNTs, suggesting that the 
TPP functionalization of multi-wall CNTs would exhibit 
the same optical properties that those observed for the 
single-wall CNTs. 

Finally, the stability of the CNT-TPP compounds is 
investigated by molecular dynamic simulations. Our re¬ 
sults show that TPP molecules in the fully covered con¬ 
figurations remain adhered on the CNT surface at 100°C, 
with adsorption distance fluctuating around 3.6 A. Fur¬ 
ther simulations suggest that the CNT-TPP compounds 
would be stable at temperatures as high as 500°C. The 
above results suggest that semiconducting CNTs fully 
covered with TPPs would exhibit excellent properties as 
a light-harvesting material. In addition, the stability of 
the porphyrins/CNT complexes is also very important in 
the field of electrocatalysis as they have shown catalytic 
activity for the O 2 reduction reactio n^^d^ , while hybrid 
metal-porphyrins/CNT systems are currently tested by 
several groups around the world as electrode materials in 
fuel cellsii. 
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